The atmospheric relevance of primary alcohols and imidogen reactions

Organic alcohols as very volatile compounds play a crucial role in the air quality of the atmosphere. So, the removal processes of such compounds are an important atmospheric challenge. The main goal of this research is to discover the atmospheric relevance of degradation paths of linear alcohols by imidogen with the aid of simulation by quantum mechanical (QM) methods. To this end, we combine broad mechanistic and kinetic results to get more accurate information and to have a deeper insight into the behavior of the designed reactions. Thus, the main and necessary reaction pathways are explored by well-behaved QM methods for complete elucidation of the studying gaseous reactions. Moreover, the potential energy surfaces as a main factor are computed for easier judging of the most probable pathways in the simulated reactions. Our attempt to find the occurrence of the considered reactions in the atmospheric conditions is completed by precisely evaluating the rate constants of all elementary reactions. All of the computed bimolecular rate constants have a positive dependency on both temperature and pressure. The kinetic results show that H-abstraction from the α carbon is dominant relative to the other sites. Finally, by the results of this study, we conclude that at moderate temperatures and pressures primary alcohols can degrade with imidogen, so they can get atmospheric relevance.

The atmospheric chemistry of alcohols is a significant global subject in the context of urban and regional air quality analysis due to their wide range of usages 1 . It has been an increased request for renewable biofuels (that are alcohol-based fuels) instead of unsustainable fossil fuels 2 . Because air pollution and climate change due to the explosive growth of the automobile industry and transportation that accelerates the consumption of fossil fuels in large quantities are global concerns [3][4][5] . On the other hand, to achieve efficient and clean combustion, the use of alternative fuels and the fuels derived from biomass which have a wide variation in their physicochemical properties carriages vast technical challenges [6][7][8][9] .
In the atmosphere and interstellar space, the most predominant non-methane organic species is methanol [10][11][12] , which is the main source of tropospheric CO 13 and formaldehyde 14 , and also has small action in the tropical HOx and ozone budgets 15,16 . Methanol is emitted into the atmosphere by several sources such as vegetation 17,18 , plant growth 19,20 , biomass burning, biofuels 21,22 , plant matter decay 23,24 , and human activities in urban and industrial regions 25 . It is better to say that in the atmosphere, methane is generated by CH 3 O 2 (peroxy radical) reactions 26,27 . The biogenic emissions of methanol have a particularly important influence on rainforest ecosystems, so they play a pivotal role in the tropospheric chemistry 28 .
Understanding the atmospheric cycle of ethanol is an important task in combustion chemistry 29 . Because it is a biogenic volatile organic compound that is highly used as a fuel for motor vehicles 30 . The same as methanol, ethanol is released into the atmosphere by different sources including vegetation and during biomass combustion. Also, various urban and industrial processes have an impact on the concentration of atmospheric ethanol 31 . The portion of ethanol that derives from biomass (bio-ethanol) is nowadays being developed as a renewable fuel that will decrease support for fossil fuels and global warming 32,33 . Moreover, there exists a beneficial profit in ethanol-utilizing as fuel for human health because of air quality 34,35 . Despite the mentioned crucial roles of ethanol that increase over time, atmospheric reports are rather sparse and the global sources and sinks for this compound are not comprehended in detail 32,36 .
Another vital alcohol is n-propanol which is considerably less toxic and less volatile than methanol 37 . Thus, it may be a desirable case for developing biofuels, but there are many discussions about that. It is well known that n-propanol serves high octane numbers and has lower corrosion than ethanol 38 . Despite these advantages, Table 1. The computed relative energies for stationary points of the CH 3 OH + NH reaction (Unit of all numbers is kcal mol −1 ). MUE 1 , MUE 2 , and MUE 3 are the mean unsinged errors for the W1BD, CBS-QB3, and M06-2X methods, respectively.   Table 2. Thermodynamic parameters for stationary points of the CH 3 OH + NH reaction (Unit of all numbers is kcal mol −1 ). A, B, and C refer to the M06-2X, CBS-QB3, and W1BD methods, respectively. www.nature.com/scientificreports/ As shown in Fig. 3, by increasing temperature, the equilibrium constants of the CR1(m) and CR2(m) computed by the CBS-QB3 method are decreased up to 800 K and after that, they are increased. For CR1(m), the obtained equilibrium constant at the M06-2X method has a similar treatment but for CR2(m), after 800 K a constant value with smooth growth is forecast. The results of the CBS-QB3 method show that the equilibrium constant of CR2(m) is 19.00, 2.34, and 1.49 times higher than CR1(m) at 300, 1000, and 3000 K.

∆E˚(A) ∆E˚(B) ∆E˚(C) ∆H˚(A) ∆H˚(B) ∆H˚(C) ∆G˚(A) ∆G˚(B) ∆G˚(C) T∆S˚(A) T∆S˚(B) T∆S˚(C)
The first path of methanol plus NH radical is highly endothermic (12.15 kcal mol −1 at the M06-2X method and 11.88 kcal mol −1 at the CBS-QB3 method) at the standard condition and has a higher energy barrier in comparison with the second pathway. In this path, NH moiety in the pre-reactive complex CR1(m) due to having a suitable orientation gets a hydrogen atom from the hydroxyl group by surmounting TS1(m) with an energy barrier of 17.61 kcal mol −1 at the CBS-QB3 level. This value is just 0.18 kcal mol −1 higher than the respective value computed by the M06-2X method. The optimized structure for TS1(m) shows that the H4 atom is located between O1 and N1 atoms (O1…H4…N1) with bond lengths of 1.345 Å and 1.120 Å, respectively. The second pathway (R 2 (m)) which has a lower barrier height is the H abstraction reaction from the methyl group. So, this is a more feasible reaction than R1(m). The difference between the energy barriers of TS1(m) and TS2(m) are 1.88 and 1.99 kcal mol −1 at the CBS-QB3 and M06-2X methods, respectively. The structure of TS2(m) accompanied by the imaginary frequency along the reaction coordinate in both methods displays that the H1 atom is transferring from the methyl group to imidogen moiety. Also, IRC calculations of TS1(m) and TS2(m) confirm the existence of the discussed pathways. The exit channels of R 1 (m) and R 2 (m) paths are met with the post-reactive complexes of CP1(m) and CP2(m), respectively. In the structure of the CP1(m), a five-membered ring structure ( ρ (RCP) = 0.0077 e bohr −3 and ∇ 2 ρ (RCP) = 0.0372 e bohr −5 ) is achieved. It contains a hydrogen bond between the atoms 1O and 4H ( ρ (LCP) = 0.0153 e bohr −3 and ∇ 2 ρ (LCP) = 0.0492 e bohr −5 ) and van der Waals bond between the atoms 1N and 2H ( ρ (LCP) = 0.0087 e bohr −3 and ∇ 2 ρ (LCP) = 0.0316 e bohr −5 ). And for the CP2(m) case, the AIM analysis demonstrates a van der Waals bond between the atoms 1N and 4H ( ρ (LCP) = 0.0298 e bohr −3 and ∇ 2 ρ (LCP) = 0.0849 e bohr −5 ).
The computed rate constants at the M06-2X (A) and CBS-QB3 (B) methods are fitted in the non-Arrhenius rate equation and depicted in Fig. 4. The following expressions are extracted in the 300-3000 K temperature range.
In Table 3, the calculated branching ratios for the methanol and NH reaction are listed. Comparing branching ratios of the P1(m) and P2(m) reveals that in the atmospheric condition, methanol degradation proceeds by hydrogen atom abstraction from the methyl group until 1000 K, and after this temperature, the OH center takes also part in removing atmospheric methanol. From energetics and kinetics points of view, it can be concluded that in the methanol plus NH reaction, the hydrogen abstraction reaction of the methyl group occurs more easily than that of the OH group. Figure 5 shows the pressure-dependent rate constants, reduced rate constants, and the ratio k ∞ /k 0 for H abstraction from the methyl group of methanol by NH. Our computed k(T,p) has positive dependence on pressure. This claim can be proved by the ratios k ∞ /k 0 and k(T,p)/k(T,1 bar) (see Fig. 5b,c) because as mentioned in the rate constant calculation section, the ratios k ∞ /k 0 and k(T,p)/k(T,1 bar) play a key role to argue the behavior of rate constants in the falloff regime. It is worth mentioning that we study the pressure effect only on the main reaction pathway of the selected reactions. For H abstraction from the methyl group of methanol by NH, the ratio k ∞ /k 0 at 600, 1500, and 3000 K, are 1.33E + 05, 5.17E + 07, and 8.07E + 08, respectively. Also, the reduced rate constant for this channel at 600 K, k(600 K,p)/k(600 K,1 bar), in p = 10 -2 , 10 -1 , 10, and 10 2 bar is 8.89E−02, 3 Tables S25-S30). In the reaction between ethanol and imidogen, three individual reaction pathways are predicted based on a hydrogen atom migration from the oxygen and alpha and beta carbons (C α and C β ) centers. The paths are summarized as follows: Each path starts with a pre-reactive collision complex. The AIM topological analysis uncovers that unlike the prereactive complexes of CH 3 OH plus NH reaction, CR1(e) and CR2(e) have ring critical points. CR1(e) has a four-membered ring structure with two van der Waals interactions. The 1N…1O interaction with the charge density of ρ (LCP) = 0.0135 e bohr −3 and ∇ 2 ρ (LCP) = 0.0540 e bohr −5 is stranger than H4…N1 interaction ( ρ (LCP) = 0.0090 e bohr −3 and ∇ 2 ρ (LCP) = 0.0344 e bohr −5 ). Also, the AIM of CR2(e) shows a ring structure with six members ( ρ (RCP) = 0.0059 e bohr −3 and ∇ 2 ρ (RCP) = 0.0256 e bohr −5 ). This complex has a hydrogen bond between the atoms 7H and 1O ( ρ (LCP) = 0.0176 e bohr −3 and ∇ 2 ρ (LCP) = 0.0532 e bohr −5 ). The CR3(e) ring is similar to CR2(e) but with five members. The hydrogen bond of CR3(e) is stronger than that of CR2(e) but the Table 3. Branching ratios of all channels in the methanol plus NH reaction. A and B refer to the M06-2X and CBS-QB3 methods, respectively. Also, S = N i=1 ki.  The graph of the computed equilibrium constants along with their fitted expressions is shown in Fig. 8. The equilibrium constants of CR1(e) and CR2(e) complexes have similar behavior compared to the equilibrium constants of CR1(m) and CR2(m) in the used computational methods. Also, the CR3(e) equilibrium constant is similar to CR2(m). The computed Gibbs free energy for CR1(e) is 5.01 kcal mol −1 which is 0.05 kcal mol −1 lower than that of CR1(m). So, CR1(e) and CR1(m) have close equilibrium constants in the 300-3000 K temperature range. Also, the free energies of CR2(e) and CR2(m) have a 0.10 kcal mol −1 difference. Thus, they have near values   Figure 5. Pressure-dependent rate constants (a), reduced rate constants (b), and the ratio of k ∞ /k 0 (c) for H abstraction from the methyl group of methanol by NH in the triplet state computed at the CBS-QB3 method. www.nature.com/scientificreports/ for equilibrium constants. In addition, this statement is correct about the comparison of CR2(e) and CR3(e) equilibrium constants (see Supplementary Table S10). The complex of the first path, CR1(e), after supplying the necessary barrier energy through the saddle point TS1(e) transforms to the corresponding post-reactive complex. The IRC calculation confirms the H atom transfer from the oxygen atom of the hydroxyl group to the nitrogen atom of the imidogen and also the existence of the minima structures of CR1(e) and CP1(e) before and after TS1(e), respectively. The obtained barrier height for the first path is 17.27 kcal mol −1 which is 0.34 kcal mol −1 lower than the barrier energy of TS1(m). Moreover, the results depicted in Figs. 1 and 6 show that the N…H and O…H bond lengths in TS1(m) and TS1(e) have negligible differences. These bonds are 0.003 and 0.002 Å longer in TS1(e), respectively. The second reaction, R 2 (e), involves H abstraction from the C α center. This reaction has a 3.77 kcal mol −1 lower energy barrier than R 1 (e). The difference in barrier energies of TS1(e) and TS2(e) may relate to the hyperconjugation phenomenon. So, the radical center on the alpha carbon can get stability by this phenomenon. Also, another reason is related to the difference in the N…H bond length in both TSs. Since the hyperconjugation phenomenon cannot take place for TS3(e), its relative energy is close to TS1(e). It should be noted that the barrier energy of TS2(e) is 5.23 kcal mol −1 lower than TS3(e). The difference is related to the abovementioned reasons and also the discrepancy of the relative energies of respective reactive complexes that CR3(e) is 0.34 kcal mol −1 more stable than CR2(e). In the  www.nature.com/scientificreports/ structure of the TS2(e), the 5H-2C covalent bond is broken up with a length of 1.346 Å, and the N-5H covalent bond is created with a length of 1.272 Å. The CP2(e) complex is the final step of this process that converts to the final product P2(e) without passing any transition state. The AIM results indicate that the post-reactive complex CP1(e) optimized structure involves a van der Waals interaction between the atoms 1N and 2H ( ρ (LCP) = 0.0077 e bohr −3 and ∇ 2 ρ (LCP) = 0.0252 e bohr −5 ), a van der Waals interaction between the atoms 4H and 1N ( ρ (LCP) = 0.0089 e bohr −3 and ∇ 2 ρ (LCP) = 0.0323), and a hydrogen bond interaction between the atoms 6H and 1O ( ρ (LCP) = 0.0135 e bohr −3 and ∇ 2 ρ (LCP) = 0.0487 e bohr −5 ), leading to form two five-membered ring structures among 1O, 6H, 1N, 4H, and 2C atoms ( ρ (RCP1) = 0.0084 e bohr −3 and ∇ 2 ρ (RCP1) = 0.0355 e bohr −5 ) and 1N, 2H, 1C, 2C, and 4H atoms ( ρ (RCP2) = 0.0067 e bohr −3 and ∇ 2 ρ (RCP2) = 0.0260 e bohr −5 ). Also, the CP2(e) has just a van der Waals interaction between the atoms 1N and 6H ( ρ (LCP) = 0.0285 e bohr −3 and ∇ 2 ρ (LCP) = 0.0819 e bohr −5 ). Also, For CP3(e), this analysis exhibits a five-membered ring structure with a van der Waals inaction between the atoms 1N and 5H ( ρ (RCP) = 0.0081 e bohr −3 and ∇ 2 ρ (RCP) = 0.0324 e bohr −5 ) and a hydrogen bond between the atoms 7H and 1O ( ρ (LCP) = 0.0165 e bohr −3 and ρ (LCP) = 0.0517 e bohr −5 ).
For the above-discussed paths, the following rate expressions are extracted from the calculated values by the TST/Eckart theory at the M06-2X (A) and CBS-QB3 (B) methods. Table 5. Thermodynamic parameters for stationary points of the C 2 H 5 OH + NH reaction (Unit of all numbers is kcal mol −1 ). A, B, and C refer to the M06-2X, CBS-QB3, and W1BD methods, respectively.
Eq.9 Eq.10 Eq.11 Eq.12 Eq.13 Eq.14 www.nature.com/scientificreports/ In Fig. 9, the graph of rate constants and associated fitted expressions are sketched. The results of rate constants for H abstraction from the OH and methyl groups of ethanol show that transferring a hydrogen atom from these groups to imidogen has near rate constants at the M06-2X method. So, roughly similar expressions 9.36 × 10 −15 T  RT cm 3 molecule −1 s −1 are extracted, respectively, for OH and CH 3 groups. But, there are small differences in the calculated rate constants by the CBS-QB3 method; i.e., the rate of the R 1 (e) path is 1.86, 2.10, and 2.00 times more than R 3 (e) at 600 K, 1500, 3000 K, respectively. These differences explicitly are seen from the rate expressions of  RT cm 3 molecule −1 s −1 (for CH 2 CH 2 OH + NH 2 formation). Also, as expected quantitatively, the ratio of the rate of R 1 (e) to R 1 (m) channels is just 1.41 (600 K), 1.27 (1500 K), and 1.23 (3000 K) at the CBS-QB3 method. Thus, in comparison with the CH 3 OH + NH reaction, there is no sensible difference in the rate of H atom transfer from the OH group of ethanol to imidogen.
(15) k 1A = 9.36 × 10 −15 T 300 3.31±0.03 exp − (11.31 ± 0.05)kcal mol −1 RT ,   RT cm 3 molecule −1 s −1 at the CBS-QB3 method are for the H abstraction reaction from the middle center of ethanol which is the main reaction channel. It should be noted that the rate of the R 2 (e) channel at 600, 1500, and 3000 K is 108.86, 8.64, and 3.84 times higher than R 1 (e) and is 204.25, 17.37, and 7.68 times higher than R 3 (e) at the CBS-QB3 method, respectively. However, the R 2 (e) channel has a 5.53, 1.94, and 1.36 times higher rate than the similar channel in methanol plus NH reaction, R 2 (m), at 600, 1500, and 3000 K temperatures, respectively.
Branching ratios of all products in the ethanol and NH reaction are collected in Table 6. According to the data in Table 6, the branching ratio of product P2 (e) in the temperature range of 300-3000 K prevails over the others. The calculated production percentage at the M06-2X method for P1(e) formation up to 1500 K is below 10% and at 3000 K is 16%. Also, a roughly similar trend is seen for the CBS-QB3 results. Instead, the branching ratio of P3(e) in both methods is different. In the M06-2X method, it has 14.7% percent of formation but in the CBS-QB3 method, the P3(e) percentage is below 10% at 3000 K.
Our compared pressure-dependent rate constants, reduced rate constants, and the ratio k ∞ /k 0 for H abstraction from the C α center of ethanol by NH are shown in Fig. 10. In Fig. 10a, the positive dependence of k(T,p) to pressure increase is seen. As mentioned above, this can illustrate more clearly by the ratio k ∞ /k 0 (see Fig. 10c). This ratio for the H abstraction from the Cα center of ethanol by NH, at 600, 1500, and 3000 K is 3.49E + 05, 1.11E + 08, and 7.24E + 09, respectively. Also, the reduced rate constants for this channel demonstrate directly this statement (see Fig. 10b). The reduced rate constant for this reaction at 600 K, k(600 K,p)/k(600 K,1 bar), when p is 10 -2 , 10 -1 , 10, and 10 2 bar is 6.28E−02, 2.90E−01, 2.21, and 3.01, respectively. The values of this ratio when p gets the same values at 1500 K are 2.20E−02, 1.58E−01, 4.87, and 1.84E + 01, and at 3000 K are 1.00E−02, 1.00E−01, 1.59, and 100, respectively. To have insight into the difference of k(T,p) between methanol and ethanol reactions, we define two ratios as (k ∞ /k 0 ) e /(k ∞ /k 0 ) m and k(T,p) e /(k(T,p) m . The prior ratio at 600, 1500, and 3000 K is 2.62, 2.14, and 8.98, respectively. And the last ratio in the pressures of 10 -2 , 10 -1 , 10, 1, and 10 2 bar and temperature of 600 K is 3.35, 3.39, 4.74, 5.81, and 6.55, respectively. These values change to 1.32, 1.44, 1.71, 1.51, and 1.47, at 1500 K and 2.02E−01, 2.37E−01, 3.09E−01, 4.93E−01, and 1.13 at 3000 K, respectively. Through the results argued in this paragraph, we understand that the rate constant of ethanol at high pressures and low temperatures Table 6. Branching ratios of all channels in the ethanol plus NH reactions. A and B refer to the M06-2X and CBS-QB3 methods, respectively. Also, S = N i=1 ki. www.nature.com/scientificreports/ is higher than methanol but at high temperatures is lower. These results show that hydrogen abstraction from the C α center of ethanol occurs easly than the methyl group of methanol when pressure increases at moderate temperatures but at high temperatures, an inverse behavior is expected.

T(K) k 1 (A)/S(A) k 1 (B)/S(B) k 2 (A)/S(A) k 2 (B)/S(B) k 3 (A)/S(A) k 3 (B)/S(B)
N-propanol plus NH reaction. The structures including geometrical parameters and PES for n-propanol plus NH reaction are represented in Figs. 11 and 12, respectively. Also, the computed relative energies and thermodynamic parameters are collected in Table 7. Moreover, other conformers and their reactions are brought out in the supplementary data (see Supplementary Figs. S3-S8, and Supplementary Tables S31-S42). The most likelihood atmospheric paths of the n-C 3 H 7 OH + NH (R(pr)) reaction are found and the main annihilation pathways are presented as follows: The same as the above-discussed paths for methanol and ethanol reactions with NH, these pathways all include conventional hydrogen atom transfer reactions.
Ring structures for prereactive complexes of n-propanol plus NH reaction are analysed by AIM theory. The same as CR1(e) and CR2(e), CR1(pr) and CR2(pr) have similar four and six-membered ring structures, respectively. The structure of CR1(pr) includes two van der Waals interactions. One of them is the 1N…1O bond ( ρ (LCP) = 0.0121 e bohr −3 and ∇ 2 ρ = 0.0458 e bohr −5 ) and the other is the C2…N1 bond ( ρ (LCP) = 0.0062 e bohr −3 and ∇ 2 ρ = 0.0208 e bohr −5 ). Also, for CR2(pr), the 8H…1O (hydrogen) bond ( ρ (LCP) = 0.0172 e bohr −3 and ∇ 2 ρ = 0.0521 e bohr −5 ) along with the N1…H5 van der Waals interaction ( ρ (LCP) = 0.0077 e bohr −3 and ∇ 2 ρ = 0.0271 e bohr −5 ) is formed before reaction. The third complex CR3(pr) has a different ring structure compared to CR3(e). The difference goes back to the van der Waals interaction that in        www.nature.com/scientificreports/ In our first suggested route R 1 (pr), the covalent 1O-9H bond with a length of 1.347 Å is ruptured and another covalent bond between the atoms 1N and 6H with a length of 1.123 Å is built. This reaction begins with CR1(pr) complex formation and after surmounting the barrier of 18.26 kcal mol −1 converts to CP1(pr). This path ends with P1(pr) product formation after the dissociation of CP1(pr) to its fragments (CH 3 CHOH and NH 2 ) without any barrier. Compare to the methanol and ethanol reactions, the hydrogen abstraction from the OH group of n-propanol is relatively difficult due to the mentioned barrier energy which is 0.65 and 1.00 kcal mol −1 higher than TS1(m) and TS1(e), respectively.

Species ∆(E + ZPE)(A) ∆E(0 K)(B) ∆E˚(A) ∆E˚(B) ∆H˚(A) ∆H˚(B) ∆G˚(A) ∆G˚(B) T∆S˚(A) T∆S˚(B)
The second, third, and fourth proposed paths are R 2 (pr), R 3 (pr), and R 4 (pr) which describe the H abstraction reaction from C α , C β , and C γ , respectively. For this, we have designed three paths that include the formation of the CR2(pr), CR3(pr), and CR4(pr) complexes due to the presence of the intermolecular interactions between reactants, the saddle points TS2(pr), TS3(pr), and TS4(pr) with the barrier energies of 10.96, 15.50, and 17.21 kcal mol −1 , and the formation of the CP2(pr), CP3(pr), and CP4(pr) product complexes. Finally, the products P2(pr), P3(pr), and P4(pr) are released after the dissociation of the bonds between the NH 2 group and the ethanol residues including CH 3 CH 2 O, CH 3 CHOH, and CH 2 CH 2 OH, respectively. The second, third, and fourth saddle points simulate the cleavage of the 3C-6H, 2C-4H, and 1C-1H covalent bonds with lengths of 1.346, 1.399, and 1.413 Å, respectively, and simultaneously the formation of the N-6H, N-4H, and N-1H covalent bonds with lengths of 1.274, 1.211, and 1202 Å, respectively.
The AIM calculations show that the CP1(pr) optimized structure has two rings that each has five members (see Fig. 11). This structure includes two van der Waals bonds between the atoms 1N and 5H ( ρ ( The rate constants of all reaction channels accompanied by their fitted expressions are depicted in Fig. 14. Also, the results of branching ratios for the n-propanol and NH reaction are tabulated in Table 8. The kinetics of the product generation via the R 1 (pr)-R 4 (pr) reaction pathways are investigated to have a suitable criterion for the conversion of reactants to each product. The fitted rate expressions for the R 1 (pr) to R 4 (pr) paths are as follows: (29)        www.nature.com/scientificreports/ Among four centers with different hydrogens from a chemical point of view, the α hydrogen has a higher rate constant. So, n-propanol destruction begins from the α carbon center. For this center, the rate constant obtained by the CBS-QB3 method is higher than that of the M06-2X method. The rate values by the expression of 1.13 × 10 −14 T     In addition, the rate constant of P3(pr) formation obtained by the CBS-QB3 method at 600, 1500, and 3000 K is 4.65, 1.10, and 0.71 times more than the generation rate of P1(pr), and is 3.32, 0.57, and 0.32 times of the rate constant of P4(pr) production, respectively. As we can see in Fig. 14, the rates of all products get close to each other at very high temperatures. Because at high temperatures, the activation energy E a = E b + mRT is being E a ≈ mRT for all reactions and the parameter m for all channels is near to each other as well. About the branching ratios, the parentage of P2(pr) generation, the same as the abovementioned reaction is high at low and moderate temperatures (see Table 8). However, the branching ratios of P1(pr) and P4(pr) products against P2(pr) and P3(pr) are negligible at low temperatures, but at high temperatures P1(pr) and P4(pr) generation are appreciable.
In Fig. 15, pressure-dependent rate constants, reduced rate constants, and the ratio k ∞ /k 0 for H abstraction from the C α center of propanol by NH are depicted. Our results uncovered that k(T,p) has a positive dependence on pressure. To express, how k(T,p) varies with pressure the ratio k ∞ /k 0 is helpful as stated before (see Fig. 15c). The calculated k ∞ /k 0 ratio for P2(pr) at 600, 1500, and 3000 K, are 3.15E + 05, 7.25E + 07, and 2.77E + 10, respectively. Also, in the pressures of 10 -2 , 10 -1 , 10, and 10 2 bar, the reduced rate constant for P2(pr) product at 600 K, k(600 K,p)/k (600 K,1 bar) Figure 15. Pressure-dependent rate constants (a), reduced rate constants (b), and the ratio k ∞ /k 0 (c) for H abstraction from the C α center of propanol by NH in the triplet state computed at the CBS-QB3 method. For the ratios (k ∞ /k 0 ) pr /(k ∞ /k 0 ) e , (k ∞ /k 0 ) pr /(k ∞ /k 0 ) m , the values 9.04 E−01 and 2.34 at 600 K and 6.54 E−01 and 1.40 at 1500 K and 3.28 and 3.43 at 3000 K are computed, respectively. Thus, in comparison with the ethanol reaction, the prior ratio has small (high) values at small and moderate (high) temperatures. But, compared to the ratio k ∞ /k 0 for methanol, the ratio k ∞ /k 0 for n-propanol has high values in the 300-3000 K temperature range. These results show that at high pressures, and low and moderate temperatures k(T,p) for propanol is higher than ethanol and methanol but at high temperatures, inverse behavior is seen. This may relate to a decrease in the reaction cross-section at high temperatures.
N-butanol plus NH reaction. Selected geometrical parameters and structures for n-C 4 H 9 OH + NH reaction are represented in Fig. 16. Also, the PES profile for this reaction is sketched in Fig. 17. The calculated relative energies along with thermodynamic parameters by the M06-2X and CBS-QB3 methods are collected in Table 9. Furthermore, we have collected the data of other conformers of n-butanol and associated reactions in the supplementary material (please see Supplementary Tables S43-S62 and Supplementary Figs. S9-S18). From a chemical point of view, butanol has different hydrogen in five positions, which may react with imidogen in atmospheric conditions. So, five different paths are designed as As the predicted paths show, by the initial association of the reactants, five different pre-reactive complexes are obtained due to different locations and orientations of imidogen around n-butanol. For starting the n-butanol plus NH reaction, five different prereactive complexes are predicted by IRC calculations. More exploration of the obtained complexes by AIM theory shows that the structure of some complexes contains not only ring critical points but also a cage critical point. The first complex CR1(bu) the same as CR1(m) has only a van der Waals interaction between the N and O atoms ( ρ (LCP) = 0.0128 e bohr −3 and ∇ 2 ρ = 0.0482 e bohr −5 ). In CR2(bu), the ring structure is the same as the rings of CR2(e) and CR2(pr) with a close value of density of ring critical points and its Laplacian ( ρ (RCP) = 0.0063 e bohr −3 and ∇ 2 ρ (RCP) = 0.0270 e bohr −5 ). The cage critical point is seen ( ρ (Cage) = 0.0056 e bohr −3 and ∇ 2 ρ (Cage) = 0.0240 e bohr −5 ) for CR3(bu) complex. This cage is formed by four intermolecular interactions. Three of that interactions are van der Waals type and one is hydrogen bond type ( ρ (LCP) = 0.0179 e bohr −3 and ∇ 2 ρ (LCP) = 0.0548 e bohr −5 ). The CR4(bu) is a seven-membered ring-like complex ( ρ (RCP) = 0.0055 e bohr −3 and ∇ 2 ρ (RCP) = 0.0222 e bohr −5 ) with two intermolecular interactions. In the CR5(bu) structure, n-butanol and imidogen moieties have two weak van der Waals interactions in a six-membered ring ( ρ (RCP) = 0.0047 e bohr −3 and ∇ 2 ρ (RCP) = 0.0160 e bohr −5 ). The equilibrium expressions of all complexes in the n-butanol plus NH reaction are as follows:   Table 9. Relative energies and thermodynamic parameters for stationary points of the n-C 4 H 9 OH + NH reaction. (Unit of all numbers is kcal mol -1 ). A and B refer to the M06-2X and CBS-QB3 methods, respectively. www.nature.com/scientificreports/ In Fig. 18, all computed equilibrium constants and their fitted expressions are displayed. The value of the equilibrium constant at the CBS-QB3 method for the CR1(bu) complex at 600, 1500, and 3000 K is 1.81E−23, 6.45E−23, and 2.57E−22 cm 3 molecule −1 , respectively. Also, the ratios of K CR1(bu) /K CR2(bu) , K CR1(bu) /K CR3(bu) , K CR1(bu) / K CR4(bu) , and K CR1(bu) /K CR5(bu) are 0.66, 0.37, 0.93, and 0.08 (600 K), 1.27, 0.82, 1.79, and 0.05 (1500 K), and 1.59, 1.08, 2.23, and 0.05 (3000 K), respectively. To compare with the previous reactions, the ratios of K CR1(bu) /K CR1(m) , K CR1(bu) /K CR1(e) , and K CR1(bu) /K CR1(pr) may be helpful. The mentioned ratios at 600 K are 1.86, 1.66, and 0.58, respectively. Our results show that the ratio of K CR1(bu) /K CR1(m) , K CR1(bu) /K CR1(e) are approximately constant and the ratio K CR1(bu) /K CR1(pr) has a small variation over 300-3000 K temperature range. For the second complexes of the discussed reactions, the value of the same ratios, K CR2(bu) /K CR2(m) , K CR2(bu) /K CR2(e) , and K CR2(bu) /K CR2(pr), at 600 K are 0.77, 1.01, and 0.68, respectively. It should be noted that the ratio K CR2(bu) /K CR2(m) remains 0.77 over the 300-3000 K temperature range. Also, other ratios change slightly at the mentioned temperature range. The ratios K CR3(bu) /K CR3(e) and K CR3(bu) /K CR3(pr) are 17.82, and 1.21 at 600 K, respectively. The same as K CR2(bu) /K CR2(m) ratio, the ratio K CR3(bu) /K CR3(pr) is constant over the 300-3000 K temperature range, but the ratio K CR3(bu) /K CR3(e) is increased to 45.91 at 1500 K and 92.64 at 3000 K. Finally, the ratio K CR4(bu) /K CR4(pr) is 0.03, 0.01, and 0.01 at 600, 1500, and 3000 K, respectively.

Species ∆(E + ZPE)(A) ∆E(0 K)(B) ∆E˚(A) ∆E˚(B) ∆H˚(A) ∆H˚(B) ∆G˚(A) ∆G˚(B) T∆S˚(A) T∆S˚(B)
In path R 1 (bu), we first study the hydrogen abstraction reaction of the OH group by NH the same as the above reactions. The energy barrier of the R 1 (bu) path is high similar to methanol, ethanol, and n-propanol reactions. In this path, the product complex CP1(bu) formation happens when a barrier of 16.86 kcal/mol is provided. The shift of the hydrogen atom from OH of n-CH 3 (CH 2 ) 3 OH to the nitrogen of the NH fragment is confirmed (42)   www.nature.com/scientificreports/ by IRC calculation. Thus, in the structure of TS1, the H-O bond is breaking up (with a length of 1.398 Å) and, the N-H bond is forming (with a length of 1.215 Å). The barrier energy of the R 1 (bu) path is 1.40 kcal mol −1 lower than the R 1 (pr) path. In the final step, the CP1(bu) complex converts to the final product P1(bu) when the interactions in the eight-membered ring structure are ruptured. The type of interactions in the mentioned ring-like structure is three van der Waals interactions and one hydrogen bond interaction. The van der Waals interactions are between the 1N and 3H atoms ( ρ (LCP) = 0.0062 e Bohr −3 and ∇ 2 ρ (LCP) = 0.0203 e bohr −5 ), the 1N and 5H atoms ( ρ (LCP) = 0.0072 e Bohr −3 and ∇ 2 ρ (LCP) = 0.0262 e bohr −5 ), and the 6H and 1N atoms ( ρ (LCP) = 0.0089 e bohr −3 and ∇ 2 ρ (LCP) = 0.0295 e bohr −5 ). And the hydrogen bond is located between the 1O and 10H atoms(ρ (LCP) = 0.0142 e bohr −3 and ∇ 2 ρ (LCP) = 0.0434 e bohr −5 ).
In the second path, the CP2(bu) product complex is obtained while alpha hydrogen of n-butanol migrates to imidogen. According to the calculated energies, this process is more feasible than the R 1 (bu) reaction due to having a lower energy barrier (13.14 kcal/mol). Also, this path has a higher energy barrier than R 2 (pr). The difference is 2.18 kcal mol -1 . In transition state 2, the length of the dissociating bond between the 2C and 4H atoms is 1.371 Å, and the length of the forming bond between the 1N and 4H atoms is 1.238 Å. The P2(bu) product is released after opening up the formed ring between the two moieties. Our results revealed that two ring critical points are situated among the 1N, 7H, 1O, 4H, and 3C atoms ( ρ (RCP1) = 0.0073 e bohr −3 and ∇ 2 ρ (RCP1) = 0.0278 e bohr −5 ) and the 1N, 5H, 2C, 5H, and 3C atoms ( ρ (RCP2) = 0.0074 e bohr −3 and ∇ 2 ρ (RCP2) = 0.0280 e bohr −5 ).
It is also possible that the reaction of gamma hydrogen observes in the gas phase. The feasible path for this hydrogen is R 4 (bu) which starts with the prereactive collision complex CR4(bu). The saddle point of this path TS4(bu) has 14.38 kcal mol −1 energy barrier which is 2.83 kcal mol −1 lower than the saddle point of TS4(pr). The transition state 4 involves cleavages of the bond between the atoms 2C and 4H (with a length of 1.371 Å) and at the same time, the bond between the atoms 1N and 4H is formed (with a length of 1.238 Å). Also, the formed post-reactive complex CP4(bu) has 4.08 kcal mol −1 lower relative energy than CP4(pr). These differences go back to the number and strength of the interactions in both complexes. For CP4(bu), AIM calculations show that it has two van der Waals bonds located between the 2C and 4H atoms ( ρ (LCP) = 0.0132 e bohr −3 and ∇ 2 ρ (LCP) = 0.0358 e bohr −5 ) and the 9H and 1N atoms ( ρ (LCP) = 0.0075 e bohr −3 and ∇ 2 ρ (LCP) = 0.0250 e bohr −5 ), causing to form a six-membered ring structure. However, after breaking the mentioned interaction, the P4(bu) product is yielded. R 5 (bu) is the last pathway that shows the simulation of the hydrogen abstraction from the methyl group with an energy barrier of 15.81 kcal mol −1 . This path after R 1 (bu) has the highest energy barrier. Also, the post reactive complex of this path after CR1(bu) has the highest relative energy. This may relate to two weak van der Waals interactions that are the reason for the formation of a six-membered ring between the NH 2 and CH 2 CH 2 CH 2 CH 2 OH moieties. Accordingly, the product of this path, P5(bu), is released in the exit channel after getting small energy (i.e., 1.48 kcal mol −1 ).
In the sequel, the reaction of n-butanol plus imidogen is evaluated kinetically. Figure 19 displays the behavior of the reaction at different temperatures. The rates of all paths are computed and collected in Supplementary  Table S16. Also, the fitted rate expressions are as follows: www.nature.com/scientificreports/ Our calculated rate constants for hydrogen atom abstraction from the OH group of n-butanol indicate that this group is less reactive at low temperatures. The fitted rate expressions at the M06-2X method,

RT
, cm 3 molecule −1 s −1 , confirm this statement due to high activation energy. It might be useful to argue about the differences in rate constants by their ratios for the same centers in different alcohols. The ratios k R1(bu) /k R1(m) , k R1(bu) /k R1(e) , and k R1(bu) /k R1(pr) at 600 K are 1.93, 1.37, and 0.64, respectively. These ratios are 1.37, 1.07, and 0.60 at 1500 K and are 1.22, 0.99, and 0.59 at 3000 K, respectively. By this comparison, we understand that there is no sensible difference among the rate constants of the OH group of simple alcohols in the reaction with NH. After the hydroxyl group, the methyl group has less reactivity. This forecast coincides with the activation energies predicted by the rate expressions of    RT cm 3 molecule −1 s −1 (CBS-QB3). These results accompanied by the above-discussed reactions lead us to conclude that the connected carbon to the OH group and also the middle carbons of the carbon chain may have higher reactivity with active species than other centers. Also, these paths have higher branching ratios (see Table 10) at moderate temperatures. Analogous to methanol, ethanol, and n-propanol, in the reaction of n-butanol plus NH reaction, the rate of Cα reaction is higher than in the other centers. This statement because of having smaller activation energies is obviously seen in the rate expressions k 2A = 5.83 × 10 −15 T  Our computed branching ratios for the n-butanol and imidogen reaction are listed in Table 10. This table shows that P2(bu) has a higher percentage of generation at lower temperatures. At moderate temperatures, P3(bu) and P4(bu) gradually become important. It should be noted that P4(bu) has a higher branching ratio than P3(bu). Above moderate temperatures and especially at high temperatures, P5(bu) also become more and more important but P1(bu) remains unimportant.
For n-butanol with NH reaction, the pressure-dependent and reduced rate constants along with the ratio k ∞ /k 0 are displayed in Fig. 20. Through Figs. 5a, 10a, 15a, and 20a, we understand that the rate of H abstraction reaction from the carbon center that is connected to the OH group is increased by pressure. For the Cα center of n-butanol in reaction with NH, k(T,p) at 600 K in the pressures of 10 -2 , 10 -1 ,1, 10, and 10 2 bar is 1.41E−17, 6.92E−17, 2.30E−16, 6.61E−16, and 5.81E−16 cm 3 molecule -1 s -1 , respectively. We reasoned this increase by the ratio k ∞ /k 0 for the above reactions. Here, we also continue the same argument. Our estimation for the k ∞ /k 0 ratio of the P2(bu) formation path at 600, 1500, and 3000 K are 3.94E + 05, 1.30E + 08, and 8.02E + 10, respectively. As Fig. 20c shows, we can say an overall remark that the ratio k ∞ /k 0 increases by carbon chain length increasing at high temperatures. This statement can explicitly be shown by the (k ∞ /k 0 ) bu /(k ∞ /k 0 ) pr , (k ∞ / Table 10. Branching ratios of all channels in the n-butanol plus NH reaction. A and B refer to the M06-2X and CBS-QB3 methods, respectively. Also, S = N i=1 ki.  It can be said that temperature increase at high pressures is slightly not suitable for n-butanol by considering the ratios k(T,p) bu /k(T,p) m , k(T,p) bu /k(T,p) e , and k(T,p) bu /k(T,p) pr (see Fig. 20a). The first ratio has values of 4.71 (at 600 K and 0.01 bar), 5.89 (at 600 K and 0.1 bar) 6.84 (at 600 K and 1 bar), 7.62 (at 600 K and 10 bar), and 7.94 (at 600 K and 100 bar). At 1500 and 3000 K this ratio has values of 1. 54  Some atmospheric events. In this section, the gas phase degradation of primary alcohols is investigated according to the real conditions of our ambient. The considering condition is the relative humidity since it plays an important role in the atmosphere. As is well proved, the main source of hydroxyl radicals in the air is water molecules. The reaction of water molecules with the produced atomic oxygen from the photolysis of ozone lead to the formation of OH radicals as follows 123,124 Also, another way for the production of hydroxyl radicals is the direct photolysis of water molecules [125][126][127] .
The formed atomic hydrogen from Eq. 59 and from other atmospheric reactions can react with ozone and generate OH radicals as well 128,129 .  Figure 20. Pressure-dependent rate constants (a), reduced rate constants (b), and the ratio k ∞ /k 0 (c) for H abstraction from the C α center of n-butanol by NH in the triplet state computed at the CBS-QB3 method. www.nature.com/scientificreports/ It should be noted that all the abovementioned reactions can be regarded as domain sources of OH production in the air [123][124][125] . So, the concentration of hydroxyl radical will change by relative humidity. For relative humidity, three states are feasible, including high, moderate, and low relative humidity. In high relative humidity, the concentration of OH radical is higher than other radical species near ~ 10 6 cm 3 molecule −151,52 . Also, theoretical 10,53-58 and experimental 11,59-68 types of research show that the rate of methanol + OH reaction has an order of ~ 10 -12 cm 3 molecule −1 s −1 at 298 K. This order for ethanol [63][64][65][66][67][68][69][70]130,131 , n-propanol 37,64,65,71,132,133 , and n-butanol 42,48,65,[134][135][136][137][138][139][140][141] plus OH reactions is almost the similar to the methanol + OH reaction. Because these rates are usually higher than the rates of other reactions of alcohols with atmospheric radicals and also due to high concentration of OH radicals, it can be concluded that the action of other radicals in comparison with OH radicals is negligible in the degradation of alcohols in an ambient having high relative humidity. But, in the relative humidity in which the concentration of OH radicals is the same as other active oxidants such as O, F, Cl, and NH species, temperature and pressure have a crucial role in the activity of atmospheric species. Thus, by considering the relative humidity of an ambient and its temperatures and pressure, and also the reactivity of atmospheric oxidants, it can refer to the following cases.
(a) In low temperatures and pressures, the activity of OH radicals is almost higher than many oxidant species (this statement is based on the reported rate constants). So, again the main sink of alcohol removal from the atmosphere is OH radicals. (b) As is well known, for many atmospheric reactions by increasing in temperature and pressure, the rate of many atmospheric reactions is increased (the inverse behavior is observed for barrierless reactions). So, in high pressures and temperatures, the rate of some reacting species such as O 142 , F 143,144 , Cl [145][146][147] , and NH is the same as OH radicals or higher. Therefore, there is a competition between F, Cl, and 3 NH radicals for the elimination of atmospheric pollutants (like alcohols). (c) In low relative humidity in which the concentration of OH radicals is negligible, other reactive species play a vital role based on their activity and concentration in the degradation of atmospheric species (like alcohols).
It has been proved that NH 3 molecules have a good concentration in the atmosphere 148 . We 108 and others [149][150][151][152][153][154] proved that the principal reaction of NH 3 is as follows: The produced amidogen (NH 2 ) molecules can undergo the following reaction 108,155 Also, amidogen radical and atomic nitrogen through the following reaction yields imidogen molecules 156 On the other hand, the photodissociation of NH 3 157-159 and HN 3 160 will help the NH generation in the atmosphere by the following dissociation reactions.
Therefore, based on Eqs. (62)(63)(64)(65)(66)(67), we understand that in high and low relative humidity, the concentration of NH radicals may be appreciable. But, as aforementioned the concentration of OH radicals in high relative humidity is very high. So, the concentration of NH is very low than hydroxyl radicals. On the other hand, our computed rate constants demonstrate that roughly at a temperature above 400 K, the reactions of NH with alcohols will take place in the atmosphere. This temperature is supplied by sun radiation mainly. Therefore, it can be understood that imidogen in dry air will act as an excellent scavenger to have a clean atmosphere from pollutants such as methanol, ethanol, n-propanol, and n-butanol. Summary. In this paper, atmospheric degradation of primary alcohols as a main atmospheric challenge was discussed because air quality directly depends on the concentration alcohols. To prove the subject, we examined the atmospheric relevance of reactions of linear organic alcohols with imidogen by using two different theoretical approaches to have comprehensive information about the chemistry of those reactions. Accordingly, two well-behaved theoretical methods, CBS-QB3 and M06-2X, based on trusted QM formalisms were used to prove the studying subject. By the used methods, all the most probable channels for removing the selected alcohols were designed and so the PES of all reactions were established. Also, the necessary energetic parameters such as thermodynamic variables and relative energies were computed for involved stationary points. Through the accomplished results, it was proven that the energy barrier of the H abstraction reaction from the connected www.nature.com/scientificreports/ carbon to the OH group (C α ) is lower than in the other centers. And, the products of this channel is more stable than the others. Thus, in atmospheric conditions, the destruction of linear organic alcohols begins through this center. To discover what temperatures and pressures in which the simulated reactions find atmospheric relevance, the temperature-and pressure-dependent rate constants were calculated for all channels of each reaction. Our computed rate constants by both QM methods revealed that the reactions of target alcohols with imidogen are meaningful at moderate temperatures and pressures. In addition, compared to methanol and ethanol, the results uncovered that the rate of long chain length alcohols and NH reactions at high pressures and moderate temperatures is high but at high temperatures, a similar rate is predicted for all alcohols. In summary, the executed procedure in this study made clear that the atmospheric relevance of the atmospheric alcohols is possible at moderate temperatures and pressures.
Computational details. For four reactions under consideration here, full geometry optimization together with harmonic vibrational frequency computations are carried out for all involved molecules such as reactants (Rs), products (Ps), and transition states (TSs) by the validated meta hybrid density functional method, M06-2X 161,162 . The used basis set for the M06-2X method is the Pople double zeta 6-31 + G(d,p) type 163 . Our recent investigations demonstrated that the structures involved in the reaction pathways and also the kinetic analysis of the H abstraction mechanism in the gas phase reactions can obtain with high accuracy by employing the M06-2X method 164,165 . After predicting reliable structures, the most popular composite method, CBS-QB3 166 , is used to calculate more precise energetic parameters. To prove the reliability of the mentioned methods 83 , the W1BD 167 method which uses the Brueckner doubles 168,169 (BD and BD(T)) methods instead of the coupled cluster 170 (CCSD and CCSD(T)) methods in the Weizmann-1 theory is applied to predict the exact energies of the stationary points. It is better to say that the W1BD is utilized only in methanol and ethanol plus NH reactions due to its high computational cost. The results revealed that the energy barriers calculated by the W1BD and CBS-QB3 methods are in satisfactory agreement with each other. According to the obtained harmonic vibrational frequencies, we find the lack of imaginary frequency in the final structures verifying that the structure is true minima. Also, the presence of only one negative eigenvalue in the Hessian (force constant) matrices confirms that the structure is a transition state (TS) one. To find the possible reaction pathways, the connectivity of saddle points to both respective per-reactive and post-reactive complexes is considered by the intrinsic reaction coordinate (IRC) calculations [171][172][173] . For IRC evaluation, the M06-2X/6-31 + G(d,p) level of theory is selected. For completing the IRC path to reach the proper minima, we optimize the geometries of the first point on the left side and the last point on the right side of the saddle point. All electronic structure calculations for extracting the geometries and energies of the components of the reaction are done by the Gaussian 09 package 174 .
The well-known theory for the topological analysis of the wave function, AIM theory, is utilized to specify the main bonding features of all structures 175 . Particularly, the electronic charge density ρ (in e bohr -3 ) and its Laplacian at critical points ∇ 2 ρ (in e bohr -5 ) analysis is used to determine the nature of the newly formed bonds. To this end, the first-order density matrix (Hessian Matrix) is computed by AIM 2000 software at the M06-2X/6-31 + G(d,p) level.
Temperature and pressure dependent rate constants calculations. The kinetics of the selected alcohols in reaction with NH is determined at different temperatures and pressures. In addition, the fitted rate expressions of the high-pressure-limit rate constants are extracted to achieve the Arrhenius parameters of all channels at various temperatures. Both abovementioned QM methods are used for calculations of rate constants over the temperature range of 300-3000 K. For pressure effect, the CBS-QBS method is utilized in the range of 1.00E−07 to 1.00E + 04 bar. The transition state theory 176 (TST) is carried out to compute the high-pressure limit rate constants of entire elementary reactions. Also, since all pathways of the above-discussed reactions involve the conventional hydrogen atom transfer, tunneling correction is important. A rather simple proposed method for computing the quantum tunneling effect is the introduced formula by Carl Eckart 177,178 . This approach is a special case of the tunneling method that is named zero-curvature tunneling (ZCT) 92,179 . The ZCT method is also referred to as ZCT-0. The ZCT-0 to establish the ground-state potential energy curve uses the Eckart function. As TST calculations predict the rate of reactions by a simple algorithm, the Eckart and ZCT approaches by using a similar manner forecast the tunneling correction. In other words, these approaches just need information on the stationary points such as reactants, transition state, and products along the minimum energy path (MEP). The Eckart method uses just the energies of the mentioned stationary points but the ZCT method in addition to the energy needs to the other parameters such as geometries, gradients, and Hessians. By these statements, one leads to conclude that both Eckart and ZCT models are the most logical choices for an acceptable predicting tunneling effect value accompanied by the TST calculations. It shall be emphasized that the produced values by these approaches are significantly more precise than the same value predicted by the Wigner method, without no extra computational attempts, when the TST theory is employed for yielding more accurate rate constant. Despite these statements, these results demonstrated that the Eckart predicted values are oftentimes higher than the values obtained by the ZCT approaches at low temperatures. And it has also been proved that tunneling is more important below a temperature of 200 K 180 . This means that the error of the Eckart method may be seen under 200 K. But, the forecasted values by the Eckart method at temperatures above room temperature are very near to the value obtained by the suggested more correct approximations 177,181,182 , including the small curvature tunneling, SCT. So, the ZCT values compared to the same SCT values are lower to some extent 181 . It seems that the origin of the error in the Eckart approach is related to the corner-cutting effects that are not included in it. However, this point is clearly noticeable that the Eckart function is often too narrow, which can compensate for the mentioned error 177  www.nature.com/scientificreports/ As we know, atmospheric reactions have different rates by variation of altitude. The change in reaction rate is mainly sensible in lower altitudes, 0.00 to 12.00 km, where for the 1.00 km change in height, about 6.49 K variation is seen in the ambient temperature. Thus, important factors such as thermochemical properties, kinetics, and so reactivity should be a function of both temperature and pressure. This conclusion leads us to investigate the influence of pressure on the reactions of linear organic alcohols to have more precise results by defining reliable conditions. Therefore, in addition to temperature, the effect of atmospheric pressure on the rate constant of the main reaction pathway of each alcohol is argued based on the chemical activation mechanism. To compute the temperature and pressure dependent rate constant, k(T,p), the strong collision approach by considering the atmospheric nitrogen (a species with high concentration) as the bath gas for producing energized molecules is chosen. The Rice-Ramsperger-Kassel-Marcus (RRKM) theory 183 by using the Ssumes program 184 is applied to predict the behavior of rate constants in the falloff regime. The used Lennard-Jones parameters for imidogen, methanol, ethanol, n-propanol, and n-butanol are 2.650 Å and 80.000 K, 3.626 Å and 481.800 K, 4.530 Å and 362.600 K, 4.55 A, and 576.680 K, and 4.688 Å and 531.3000 K, respectively 185 . The energetic parameters of the CBS-QB3 method and the Lennard-Jones parameters are used to generate the input of RRKM calculations. It should be noted that k(T,p) is calculated in the pressures of 1.00E−7 to 1.00E + 4 bar over the 300-3000 K temperature range. For H abstraction channels, k(T,p) graphs of methanol, ethanol, n-propanol, and n-butanol reactions with NH are sketched in Figures, 5,10, 15, and 20, respectively. Also, k(T,p) values of the mentioned reactions are listed in Supplementary Table S17-S20, respectively. The simple form of the pressure-dependent rate constant is as follows 164,183 : As we can see from the simple form of the equation defined for the pressure-dependent rate constant, the ratio of k ∞ /k 0 is a key term to forecast the falloff regime behavior of the rate constant. Also, defining the reduced rate constant as the ratio of k(T,p) into k(T,1 bar) gives a piece of clear information about the effect of pressure on the reaction of each channel. Therefore, through both ratios of k ∞ /k 0 and k(T,p)/k(T,1 bar), we consider the pressure effect on the gas-phase reactions of methanol, ethanol, n-propanol, and n-butanol with NH.
A standard software, Gpop program 186 , is implemented for computing all high-pressure limit rate constants.

Data availability
All data generated or analysed during this study are included in this published article [and its supplementary information files].